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Abstract 

In this paper we study the motion of two particles diffusing on low- 
dimensional discrete structures in presence of a hard-core repulsive in- 
teraction. We show that the problem can be mapped in two decoupled 
problems of single particles diffusing on different graphs by a transforma- 
tion we call diffusion graph transform. This technique is applied to study 
two specific cases: the narrow comb and the ladder lattice. We focus on 
the determination of the long time probabilities for the contact between 
particles and their reciprocal crossing. We also obtain the mean square 
dispersion of the particles in the case of the narrow comb lattice. The 
case of a sticking potential and of 'vicious' particles are discussed. 
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1 Introduction 



The diffusion of interacting particles has recently become the object of a renewed 
interest due to the great variety of its applications [1-14]. The first systematic 
approaches to the problem were introduced in mathematical literature by Harris 
[0 and then by Fisher in Q , where the simultaneous motion of p random walkers 
was studied and multiple occupation of a single site was forbidden by the hard 
core repulsion. Since then, different kinds of interactions have been considered, 
like the so called 'n-friendly' walkers (two of the walkers can move together for 
up to n lattice sites) and the 'vicious' walkers (walkers kill each other when they 
meet and the diffusive process stops) ^ [| [| . 

In general, the case of two diffusing interacting particles can be regarded as 
the motion of one particle moving on a randomly evolving environment. This 
gives rise to highly non linear behaviors and implies a strong influence of the 
underlying geometry on the process. This has already been evidenced in the 
above mentioned works, even if only standard d-dimensional lattices have been 
considered. 

Interactions play a dramatic role in low dimensional systems, where multi- 
ple collisions between particles are so frequent that they completely modify the 
usual random walks behavior. Extensive studies have been carried on in the 
simplest case of 1-dimensional lattices where anomalous behaviors have been 
evidenced [gj. Now since relevant applications of diffusion with interactions 
concern biological systems, complex networks and abstract spaces which typi- 
cally exhibit low dimensionality but non standard geometry, the usual modeling 
with regular lattices is not appropriate. These structures are naturally repre- 
sented by graphs and their study requires new mathematical approaches based 
on graph theory. 

In this paper, we focus on the diffusion process of two particles in presence of 
a contact repulsive interaction on low-dimension structures of non conventional 
geometry called the 'Narrow Comb' and the 'Ladder' lattice. These graphs 
reproduce the geometrical feature of some simple polymers, like polymeric liquid 
crystals, and decorated linear networks. Though keeping the same large scale 
structure of a linear chain, they allow particle crossing in spite of the presence 
of a repulsive interaction. 

The problem of two diffusing particles on these graphs is solved introducing 
a new technique we call the diffusion graph transform (DGT) fl5|| . The original 
problem can be decoupled into the motion of two independent random walkers 
moving on two graphs, representing the coordinate of the center of mass and 
the relative distance between the two particles respectively. Then the problem 
can be analytically solved using discrete time random walks methods on non- 
translation invariant lattices |l6|, |17fl . Moreover, modifications of the potential 
(sticking potential or vicious walkers potential) can be easily taken into account. 
The paper is organized as follows. In section 2 we present DGT in the simple 
case of the regular 1-dimensional chain. In section 3 we study and solve the case 
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of the Narrow Comb while on section 4 we consider the Ladder graph. Finally 
in section 5 we discuss our results. 



2 The diffusion graphs 

As a first example of application of the diffusion graph technique, we present 
the simple and well studied case of the linear chain. Let us consider a 1-d lat- 
tice with two classical distinguishable particles with coordinates x\ and x 2 . At 
t = the particles start from two adjacent sites x\ = and x 2 = 1. At each 
discrete time step they can jump to one of their nearest-neighbors with the same 
probability p — \ if these sites are unoccupied. When the particles are on two 
adjacent sites, they are forced to move apart and to occupy their only empty 
nearest neighbor: these prescriptions model the hard-core repulsive interaction. 
One of the main features on the linear chain is that the two particles cannot 
cross each other and therefore x\ < x 2 - 

The key step in the diffusion graph technique is the introduction of the coordi- 
nate of the center of mass 



x 2 +xi 1 
2 2 



(1) 



and the value of the relative distance 



xi- x\ 1 

dr = —5— ^ (2) 

where the additional term — ^ defines the initial conditions c m = d r = 0. 
At each discrete time step, the value of c m can change to c m + 1 and c m — 1, 
corresponding to both particles moving to the right or to the left, or can remain 
unchanged if particles move in opposite directions. The probability of these 
moves can be easily calculated from the original process rules and they can be 
shown to be respectively p = 1/4 and p = 1/2. An analogous derivation can 
be obtained for the value of d r , leading to the same conclusions, with the only 
additional prescription that d r cannot assume negative values. This suggest to 
represent the value of c m and d r as the positions of two random walkers moving 
on two new lattices we shall call the diffusion graphs. 

In the case of c m the diffusion graph will be a linear chain with jumping proba- 
bilities p =1/4 and waiting time probabilities in each sitep = 1/2 as represented 
in Fig.l a). The diffusion graph of d r , represented in Fig.l b), is a half linear 
chain with waiting time probabilities p = 1/2 in every site, except in the origin 
O, since this point represents the contact between the particles (d r = 0). Once 
the diffusion graphs are defined, the original problem is reduced to the study of 
a single random walker on the diffusion graphs of c m and d r respectively. 
In each single walk the time evolutions of c m and d r are correlated in such a 
way that, when one of these coordinates changes, the other is left unchanged. 
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However this very particular kind of correlation has no influence on the average 
values since it is completely taken into account by the waiting time probabilities 
on each diffusion graphs. The only residual correlation appears in a modified 
waiting time probability when the two particles have a contact, since in that 
case d r must increase of one unit, with a corresponding waiting time probability 
p = 1 for c m instead of the usual p = 1/2. However it can be proven that this 
change in the waiting probability does not affect the leading asymptotic behav- 
iors for t — > oo Jlj|. 

By standard random walks technique on the diffusion graphs, one obtains the 
relevant quantities for two interacting particles on a chain. The probability of a 
contact between particles is mapped into the probability of returning to point 
O for a random walker on the d r diffusion graph and follows the asymptotic 
law: Po(t) ~ l^/x^ 1 ^ 2 as t — > oo. The coordinate of each particle and of the 
center of mass grows as t 1 / 2 . ||]. We can also consider a sticking probability 
after a contact simply by introducing a waiting time probability in point O of 
the d r diffusion graph. Finally, the problem of 'vicious particles', killing each 
other when occupying nearest neighbors sites, can be described by considering 
only the probabilities of reaching O for the first time on the diffusion graphs. 

3 Two interacting particles on the narrow comb 
lattice 

The possibility of particles crossing is introduced on the linear chain by adding 
fingers of unit length at each site. The corresponding geometrical structure is 
the narrow comb lattice (NCL) represented in Fig. 2 a. 

Let us consider two particles starting at t = from the adjacent sites x\ = and 
x-2 = 1 on the backbone of the comb. The particles move on the NCL according 
to the following rules: 

i) each particle can jump from site i to one of its nearest neighbors sites with 
probability l/zi, Zi being the number of nearest neighbors of site i. This means 
that when a particle is on one of the fingers, it is forced to jump to the backbone 
with probability equal to one. 

ii) When the two particles occupy two adjacent sites on the backbone (contact), 
the particles are forced to move in one of the following ways: 

a) the two particles jump on the backbone following opposite directions 

b) both particles jump on the fingers corresponding to their positions 

c) only one of the particles jumps on the finger while the other moves on the 
backbone. This can be done in four different ways. Two of them lead to what 
we shall call a 'tower configuration', where the particles occupy a site on the 
backbone and the site of the corresponding finger. 

Hi) the only allowed evolution of a tower configuration consists in the lower 
particle taking a further step on the backbone, while the particle on the finger 
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drops down to the backbone. This means that particles have a probability 1/2 
to cross each other, depending on the lower particle taking the step to the left 
or to the right. 

Now we can build the diffusion graphs for the center of mass c m and the relative 
distance d r , with c m and d r defined in (Q) and (Q), where x\ and x 2 are the 
projections on the backbone of the positions of particles 1 and 2 respectively. 
When one of the particles is on a finger, c m and d r have non integer values and 
d r has negative values when x\ > x 2 - The diffusion graphs of c rn and d r are 
themselves narrow comb lattices, decorated with loops and triangles as in Fig. 
3 and Fig. 4. The position of c m on the NCL can be easily deduced from the 
position of a random walker on the corresponding diffusion graph: when the 
random walker occupies a site on the backbone, its position corresponds to the 
value of c m ', when it stays on a finger, the value of c m is that of the position of 
that finger and both particles occupy sites on the fingers of the NCL. Finally, 
when the walker is on the upper vertex of one of the triangles, c m has a semi- 
integer value corresponding to the projection of the vertex on the backbone. 
In this situation one of the particles is on the backbone and the other on the 
fingers. Loops are introduced on the diffusion graph of the NCL to represent 
waiting time probabilities for c rn . 

The relative distance graph can be built in a similar way, taking contact into 
account. Point O represents the tower configuration which can evolve on point 
A (contact with x\ < x^) or in point A' (contact with x\ > x 2 ) and all the other 
points will be indicated with the value of their projection on the backbone (1/2, 
1, 3/2, 2... and —1/2, —1, —3/2, —2..., depending on the walker being on the 
left or the right side of the graph with respect to point O). 
The diffusion graph with the corresponding jumping probabilities is represented 
in Fig. 4. The jumping probabilities not explicitly indicated are equal to the 
corresponding ones of the c m diffusion graph. Since points A and A' represent 
contacts between particles, we can have particle crossing if A — > A' or A' — > A 
i.e. crossing can only be the consequence of a contact. Now we can study the 
motion of the two interacting particles analyzing random walks on the diffusion 
graphs. 

The first quantity we calculate is the probability of a contact between parti- 
cles without crossing. This corresponds to the probability that a single walker 
starting from point A of the diffusion graph of d r returns to point A after t 
steps without ever reaching point A'. We call this quantity P^ A \A,A;t) and 
we define F^ A ) {A, A\ t) as the probability that a walker starting from point A 
returns to this same point for the first time after t steps without ever reaching 
point A'. We have: 



ipW(l,l;t-2)- 
6 









-;t- 2 ) 
2' J. 



(3) 
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where the first term refers to the two steps A — > O — ► A, the second to the 
probability of reaching the finger corresponding to point A and then going back 
to the backbone. p( n '(i,j;t) is the probability of a t steps walk from point i to 
point j never reaching point n. Let us introduce the generating function of a 
generic P( n >(i,j;t) as: 



(4) 



t=o 



Using the relation P {n \i, i; A) = 1/[1 - i; A)][|§, the calculation of one 

of the P( n '(i,i; A) is reduced to that of the corresponding F^ n '{i,i] A) which is, 
in general, a much easier task. In the case of F^ A ) (A, A; t) we obtain: 



F {A '\A,A-\) 



A 2 A 2 A 2 



ipW(l,l;A)+ 



ip<*> ( i, i; a") + ipw ( i, I x) + f i, i; A 



2' 2' 



(5) 



The generating functions of P^(l, 1; A) and P^(l/2, 1/2; A) can in turn be 
expressed as functions of the corresponding first time generating functions: 



F(^)(1,1;A) = 7 |a 2 — 



27 1-A 



2' 2 



27 1 -x 



(6) 
(7) 



where: 



^(i,i;A) = ^ + Ha + ^ 



-P«(2,2;A)- 



9 V 2 / 3 V 2 / 3 V 2 2 /. 
Moreover since P^ n \i,j; A) = F^ n '(i,j; X)P( n '(j,j; A) we also have: 



:,1;A 



1-| 3 



P^(1,1;A) 



V 2' / 1-x 9 V 2 2 / 



(8) 



(9) 



(10) 



Using the translation invariancejn the backbone direction of the^comb, we get 
the equalities: P«(2,2;A) = PW(1,1;A); P«(3/2,2;A) = P^(1/2,1;A); 
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P«(2,3/2;A) = P( A )(1,1/2;A) and pW (3/2,3/2; A) = P^(l/2, 1/2; A). We 
obtain from (0) an equation for x: 



A 2 2 A 2 



3 + 3A + 18(l-x) 



(11) 



9(3- A)(l -x) -2A 2 
Now we can substitute the value of x and obtain the hnal expression: 

(A A;X) = — + — (l+ f + 3A + 18(1 --M (12) 
v ' 7 6 6 V 9(3-A)(l-x)-2A 2 J y ' 

The asymptotic behavior of P^ A \A,A;t) can be derived from that of the cor- 
responding generating function considering the limit A — > 1 and then applying 
Tauberian theorems fllffl: 



P {A '\A,A;t) ~ 12y^t~4 



i^oo (13) 



This asymptotic law shows a dramatic difference with respect to the case of 
the simple linear chain where particle crossing is forbidden and the probability 
of a contact decays as t^ 1 / 2 . Notice that the deep change in the asymptotic 
behavior, originated by the new geometry, has no analogous if we consider the 
motion of a single particle, since in this case the linear chain and the NCL have 
the same asymptotic diffusion laws. 

In a similar way, we can calculate the probability of particles crossing. This is 
the probability P(A,A';t) that a random walker starting from point A of the 
diffusion graph of d r reaches point A' after a t steps walk. Since 

P{A, A'; A) = F(A, A'; X)P(A', A'; A) (14) 

we find after some calculations: 

P(A,A';t)~~^t-i t^oo (15) 

A key quantity is the average number of contacts without particles crossing. 
This is given by the average number of returns to point A after t steps on the 
d r diffusion graph without reaching point A': 

M {A '\A,A]t) ~ const (16) 

Finally, the mean number of particles crossing follows the asymptotic law 

M(A,A',t) ~ \\J^£ t^oo (17) 

For a better description of the diffusion properties of the two particles on the 
NCL, we study the mean square dispersions Ax 2 and Ax 2 , defined as 

Ax 2 = Axj =< x\ > - < xi > 2 =< x\ > - < x 2 > 2 (18) 
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In H the mean square dispersions of two hard core interacting particles on a one 
dimensional lattice was studied. In this case, particle crossing was forbidden by 
the potential and by the geometry of the system and it was found that 

Ax\ = Ax\ = (l-^J ct (19) 

c being the diffusion constant, while for a single particle moving freely on a one 
dimensional lattice we have 

Ax\=ct (20) 

From the previous results it follows that the interaction inhibits through the 
factor (1 — 1/tt) the spreading. In the NCL case, from ([j]) and (0) it follows 
that: 

< xi >=< c m - d r > (21) 

< x 2 >=< c m + d r > +1 (22) 

The expressions of < c m > and < d r >, due to the diffusion graph transform, 
can be calculated as the mean displacement < d > of a single walker moving 
on the correspondent diffusion graph, d being the coordinate of the projection 
of the position of the walker on the backbone of the diffusion graphs. Using the 
equation 

n^noAt) (23) 

and following the steps described in |17j we find < c m >=< d r >= 0: this also 
immediately follows from the fact that the diffusion graphs of c rn and d r are 
symmetric with respect to point O so that < c m > and < d r > must be 0. 
Equations (^l]) and ( p2| ) become 

< x-y >=< c m - d r >= (24) 
< x 2 >=< Cm + d r > +1 = 1 (25) 

and we have also 

Aa^ = Axj =<c 2 rn > + <d 2 r > (26) 
By straightforward calculations we find 

< <& >~ ^ f - oo (27) 

< d 2 r >~ -t t^oo (28) 



so that 
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Axj = ~ 2* 1 ~* 00 ( 29 ) 
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For a single particle diffusing on the NCL we have 

Ax 2 ~ -t t — ► cxd (30) 

where x is defined as the projection of the particle position on the backbone of 
the lattice (for a single walker on a linear chain Ax 2 ~ t t — > oo). We can 
therefore conclude that the interaction between particles does not modify the 
spreading as a consequence of the possibility of particle crossing. This result is 
confirmed by equations (|l^) and Jl5|), stating that a contact with crossing is 
much more frequent than a contact without crossing. 



4 Two interacting particles on ladder lattices 

The ladder lattice exhibits new features with respect to the previous cases due 
to the presence of loops. The ladder graph is made of two linear chains whose 
corresponding points are connected by links as in Fig. 2 b. The two particles 
start from the same site and then move on the ladder with jumping probabilities 
equal to 1/3. The relative distance is defined as d r = N/2, N being the chem- 
ical distance i.e. the number of links of the shortest path connecting the two 
particles. The d r diffusion graph is represented in Fig. 5; the value of d r is given 
by the projection on the lower chain of the position of a random walker. When 
the walker is on the lower chain of the d r diffusion graph, particles on the real 
ladder lattice are on the same chain, while when the walker jumps on the upper 
chain of the d r diffusion graph it means that particles move on different chains. 
Point O represents the contact between particles and the link between A and A' 
represents the possibility of particles crossing each other without direct contact, 
since it connects points with d r > and d r < without passing through point 
O. This is the direct consequence of the presence of loops and it represents a 
fundamental technical difference with respect to the previous case. 
The solution of the random walks problem on the diffusion graph can be sim- 
plified by the following considerations. As already shown in the diffusion 
problem on a ladder-like lattice can be mapped into that of a linear chain. This 
follows from the fact that we are interested only in the walker displacement in 
the unbounded direction, and not in its position on the upper or on the lower 
chain. Then the ladder can be considered as a linear chain where we define ad- 
ditional staying probabilities corresponding to the steps spent for jumping from 
one chain to the other. In the present case this simplification can be applied 
to the two half-linear chains originating from the central loop. Notice that the 
central loop cannot be reduced to a couple of points of the linear chain and 
must survive our transformation since it breaks the translational symmetry of 
the network. This detail is crucial for the calculation of crossing without contact 
probability. 

The final graph obtained by this analysis is represented in Fig. 6. It reproduces 
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the main features of the d r diffusion graph for the probability of particles cross- 
ing with or without contact. From calculations analogous to the previous case 
we find that the probability of contact without crossing is: 

pM(A,0;t)~^Jlt-l t^oo (31) 

and shows the same time dependence we found in the NCL case. The probability 
of crossing without contact is: 

P«»(A,A';t)~^Jlt-* t^oo (32) 

za V 7T 

while the generic probability of crossing decays as: 

3 i 

P(A, A'; t) i=t~* t -> oo (33) 

v 2n 

The description of the steps leading to (|3l|), (p^), (|3^ ) can be found in Appendix 
A. 

We see that the probability of a contact without crossing on the ladder lattice 
follows the same law we found in the NCL case and this is also true for the 
crossing probability. On the ladder lattice however, a probability of crossing 
without contact can be defined. This particular feature of the ladder lattice can 
be extremely useful when the two processes representing contact and crossing 
are used to represent different physical situations. 

5 Summary and Discussion 

The diffusing graph transform introduced in this paper is a powerful technique 
which, combined with the universality properties of the long time asymptotic 
behavior of random walks return probabilities, provide a new tool for the anal- 
ysis of interacting diffusing particles on discrete structures. The technique can 
be applied to decouple and analytically solve the problem on tree-like and quasi 
1-dimensional structures also in presence of loops and crossing probabilities for 
the particles. This has been done for the case of two interacting particles on the 
Narrow Comb and on the Ladder lattice. Our solution shows that the asymp- 
totic behavior of contact probabilities is sensitive to the local geometrical details 
of the underlying discrete structure. 

Appendix A 

The explicit calculation of ([3j]), (f32|), ( |33|) is based on two properties of the 
generating functions |L6]]. The first states that the generating function of the 
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probability P^(i,j;X) of going from point i to point j (i ^ j), without ever 
touching point k, is related to the corresponding first time arrival generating 
function by the relation 



pM(i,j;\) = FW(i,j;\).pW(jJ;\) 



(A.l) 



The second property is the well known relation between P^ (i, i; A) and F^ k > (i, i; A) 
we already used in the preceding sections; 



P«(M;A) 



(A.2) 



1- F^(i,i;X) 

Let us start from the calculation of P^ A \A, 0; A): from (A.l) it follows that: 
P {A '\A, 0; A) = F^ A '\A, 0; A) ■ P^ A '\0, 0; A) (A.3) 



since: 



and 



F^ A '\A,0;X) = ^P^ A '\A,A;X) 



F (A,) (0,0;A) 



1 -F^')(o,0;A) 1- %P(°> A 'HA,A;X) 



(A.4) 



(A.5) 



so that the calculation of P I ^ A '\A, 0; A) is reduced to that of P [ °- A '\A, A;X). 
To obtain the expression of p(°' A ) (A, A; A) we calculate the corresponding first 
time arrival generating function: 



F^ A '\A,A;\) = ^+(^) P^(1,1;A) 



2A 



(A.6) 



From the translational invariance of the lattice, we have that P^(1,1;A) = 
p(0,A )^ ^. ^ Using the preceding relation together with (A.2) we obtain: 



p(0,A')( AiA . X) . 



5 A 



A 2 10A 



1 + 



s.\- 



(A.7) 



In the limit A — > 1 — e and e^Owe have: 



P { °- A '\A,A-1- e) 



■x 1-3W- 



and: 



P^(A0;l-e)^-|^ 



(A.8) 
(A.9) 
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Applying Tauberian theorems we obtain @. For the determination of P (0) (A, A'; A), 
we follow the same technique writing this quantity as a function of p(°' A > [A, A; A) , 
taking the limit A — > 1 — e(e — > 0) and using Tauberian theorems. In conclusion 
we have: 

F {0 \A,A';\) = —?(°> A '\A,A;\) (A.10) 
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(A', A'; A) = 1 - ^ +(—) P { °- A ' ] (A, A; A) 



(A.ll) 



P(°^')(A,A; A) 
and finally using (|A.2| ) we have 

P (0) (A, A'; A) = F (0) (A, A'; A) • P (0) (A', A'; A) (A.12) 

For P(^4, A'; A) the corresponding steps can be resumed as: 

P(A, A'; A) = F(A, A'; X)P(A', A'; A) (A.13) 

F(A', A'; A) = 1 - ^ + 

A 2 



12- \ 2 P(°> A ')(A,A; A) 



(lp^')(4,A;A) + ^°^)(Ai;A) + l 



(A.14) 

F(A, A 1 ; A) = ^°f(^A;A) /2A A 
V ^ 12- A 2 F(o^')(A,A;A) V 3 / V 7 
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Figure captions 



Fig.l 

The diffusion graphs of a linear chain: a) c m diffusion graph and b) d r diffusion 

graph 

Fig.2 

a) 'Truncated comb lattice' b) 'Ladder lattice' 
Fig.3 

c m diffusion graph and jumping probabilities for the truncated comb lattice 
Fig.4 

d r diffusion graph and jumping probabilities for the 'Truncated comb lattice' 
Fig.5 

d r diffusion graph and jumping probabilities for the 'Ladder lattice' 
Fig.6 

Semplified d r diffusion graph for the Ladder lattice 
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